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We theoretically consider the amplification of THz radiation in a superlattice Bloch oscillator. The 
main dilemma in the realization of THz Bloch oscillator is finding operational conditions which allow 
simultaneously to achieve gain at THz frequencies and to avoid destructive space-charge instabilities. 
A possible solution to this dilemma is the extended Limited Space-Charge Accumulation scheme of 
Kroemer (H. Kroemer, cond-mat/0009311). Within the semiclassical miniband transport approach 
we extend its range of applicability by considering a difference in the relaxation times for electron 
velocity and electron energy. The kinetics of electrons and fields establishing a stationary signal in 
the oscillator is also discussed. 
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I. INTRODUCTION 

Terahertz radiation (0.3 — 10 THz) has enormous 
promising applications in very different areas of science 
and technology such as space astronomy, wideband com- 
munications and biosecurity, to name a few (for recent re- 
views, see pj). One of the main challenges is to construct 
a coherent miniature source of THz radiation that can 
operate at room temperature. Currently many groups 
worldwide are working in this direction. Along with 
many experimental and technological problems there are 
still several fundamental problems to be solved. In par- 
ticular, the traditional lasing scheme is based on a popu- 
lation inversion between different energy levels. A great 
recent achievement was the development of quantum cas- 
cade lasers that can operate at THz frequencies employ- 
ing quantum transitions between the energy levels in 
multiple quantum well heterostructures Q . These quan- 
tum nanodevices require a low temperature to achieve 
a significant population difference. Continuous improve- 
ments in the design of THz quantum cascade lasers allow 
to increase the temperature of operation above 100 K 
However, it would be a quite difficult, if ever pos- 
sible, to reach population inversion at room tempera- 
ture. Really, the spacing of the energy levels, which is 
necessary to emit radiation at frequencies of the order 
of 1 THz, is comparable with the room temperature in 
proper units. This simple observation makes it attractive 
to study other suggestions for nano-devices that do not 
require population inversion for the generation of THz 
radiation. 

The Bloch oscillator is an inversionless THz laser that 
is based on Bloch oscillations of miniband electrons in 
a dc-biased semiconductor superlattice. The continuous 
operation of a Bloch oscillator has still not been demon- 



strated in any experiments. The development of a super- 
lattice Bloch oscillator has been initiated by suggestion of 
Esaki-Tsu Q and theoretical analysis of Ktitorov, Simin 
and Sindalovskii The main obstacle in the experi- 
mental realization is the formation of high-field electric 
domains in superlattice. 

Nowadays the problem of THz Bloch oscillator again 
attracts much attention (for review, see [E0.0)- There 
are several interesting suggestions for modifications of 
the original scheme of Bloch oscillator, which in principle 
should allow to reach THz gain in superlattice without 
formation of destructive electric domains. 

In this report, we first briefly review the static electric 
properties of superlattices and their influence on both 
small-signal gain and electric instability in section [H] In 
the original part of this report we develop the idea of a 
large-amplitude THz Bloch oscillator suggested by Kroe- 
mer [9|. We consider a dc-biased superlattice subjected 
to a monochromatic ac probe field, so that the total elec- 
tric field acting on the miniband electrons is 



E = E, 



dc 



Ear COS Ujt. 
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Note that in a real device the probe field with the am- 
plitude E ac should be a mode of resonator tuned to a 
desirable THz frequency. It is well known that while a 
superlattice device should allow strong gain for a small ac 
field (E ac — > 0) in the wide range of frequencies from zero 
up to several THz, simultaneously destructive electric in- 
stabilities would arise inside superlattice (see section ITT)l . 
As Kroemer has demonstrated for particular values of 
superlattice parameters, the ac field with a large enough 
E ac can suppress domains but still preserve significant 
THz gain p. In ref. the single scattering constant 
approximation was used. However, a more realistic ap- 
proximation allow two different relaxation constants, 7„ 
for the electron velocity and 7 e for the electron energy. 

Using the technique of superlattice balance equations 
( section ITTTT) . we re-examine the THz gain and the crite- 
rion of electric stability for the case of large E ac . Our 
main aim is to find how robust the Kroemer scheme of 
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Bloch oscillator is against the effect of different relax- 
ation constants (section ^VJ. Here our main finding is 
that this effect does not dramatically change the results 
obtained within the single scattering constant approxi- 
mation, if ratio j £ /lv > 0.1. Moreover, for ~f £ /lv > 0.5 
the differences from the results obtained within the single 
relaxation time approximation are insignificant. 

In section we consider the characteristic time scales 
for the development of space-charge instabilities and for 
the growth of the ac field due to high-frequency gain 
in superlattice. We come to the conclusion that it re- 
ally is possible to suppress a destructive accumulation of 
charges in the case of THz oscillations with large am- 
plitudes. However, our estimates still demonstrate that 
it is very difficult to switch the device in the regime 
of large amplitudes before electric domains would be 
formed. Discussion devoted to this and other remain- 
ing problems in the realization of THz Bloch oscillator is 
presented in the final section. 



II. STATIC ELECTRIC CHARACTERISTIC 
AND SMALL-SIGNAL GAIN 

Let us first review the static electric properties of su- 
perlattices using a single relaxation time. Nonlinear elec- 
tron transport in a superlattice with applied dc bias E^ c 
is a well-studied problem. The dc current through the 
superlattice depends on the dc bias as Q 



Iet — I* 



peak 
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Here E cr — h/(edr) is the Esaki-Tsu critical field, 
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is the peak current corresponding to Ed c = E cr , r is the 
scattering time, n is the density of carriers, vq — Ad/ (2h) 
is the maximal electron velocity in the first miniband, A 
is the miniband width, d is the period of superlattice, 
T is the temperature, A is the cross sectional area of 
superlattice and I\ (x) and 1$ (x) are the modified Bessel 
functions. 

The dependence of Iet on Ed c is shown in Fig. ^ 
For Edc. > E cr the current-field characteristic demon- 
strates negative differential conductance (NDC). Note 
that Edc/ E cr — war with lob = edEd c /Ji being the Bloch 
frequency. Therefore the condition of static NDC is also 
ljbt > I- For typical semiconductor superlattices and 
for the applied dc bias ~ 1 kV/cm, the Bloch frequencies 
belong to THz range Q. 

In the following it is useful to consider the Esaki-Tsu 
current as a function of voltage drop over one super- 
lattice period 



lET(eEd c d) = -/pcak" 



2eE dc d/T 
f {eE d cd/Tf 



(3) 



with r = Ti/t. 

We turn to the consideration of electron transport un- 
der the action of combination of dc and ac fields. The 
time-dependent current induced in superlattice by the 
field can be represented as 

oo 

i it) = io + E[ j r cos cos M + c sin smM]. (4) 

h=\ 

Absorption of the probe ac field is proportional to the 
first Fourier component of the time-dependent current 
Ii° s . Negative absorption or gain corresponds to If° s < 
0. 

We start the discussion of small-signal gain (E ac — > 0) 
with the case of quasistatic interaction lot <C 1. In this 
case, it is easy to show that If os is proportional to the 
slope of static VI characteristic Therefore, a choice 
of working point in the NDC portion of superlattice VI 
characteristic results in small-signal gain for a quasistatic 
ac probe field. Obviously, in the quasistatic case the gain 
does not depend on the frequency. For typical superlat- 
tices at room temperature r ~ 100 fs. Therefore the 
quasistatic approximation is valid for microwave fields. 

For a higher frequency (lot > 1), it is still possi- 
ble to have small-signal gain, but the magnitude of the 
gain becomes frequency-dependent. Small-signal absorp- 
tion is also often defined using the frequency-dependent 
complex conductivity (j(lo): If os cos(wi) + Jf m sin(ajt) = 
5ie {a(Lo)E ac e~ luJt A] . It is easy to see that absorption is 
proportional to 3?e o(lo). Complex high-frequency con- 
ductivity has been first calculated by Ktitorov et al. 
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where a sta t = Iet / {AEd c ) is the static conductivity of 
the superlattice. Real and imaginary parts of the conduc- 
tivity 10 are shown in Fig.|2J Importantly, for Ed c 3> E cr 
gain exists up to frequencies around the Bloch frequency, 
with even a resonance near lo b ■ The physical mechanisms 
of this high-frequency gain has been studied both in the 
semiclassical picture with the help of electron bunches 
|lO| and in the quantum mechanical Wannier-Stark pic- 
ture [13 • 

Alternatively, the small-signal absorption can be cal- 
culated by taking the so-called quantum derivative of the 
current-field characteristic 0, 0] (Fig. ^) 



rcos 
1 l 



lET(eE dc d + Tiu) - I E T{eE dc d - tuo) 
2~Klo 



eE ac d. (6) 



For lot the difference quotient in the quantum 
derivative becomes the usual derivative, and the re- 
sult that small frequency absorption is proportional to 
<IIet / dE dc is rediscovered. One immediate consequence 
of the quantum derivative is that in order to get high- 
frequency gain we must have static NDC (see Fig. [Q. 
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FIG. 1: (Color online) Esaki-Tsu current-field characteristic. 
The figure also demonstrates the quantum derivative for the 
frequency u) = 3r _1 taken at the operation point Edc = 5E cr 
(blue). Absorption at frequency u) is proportional to the dif- 
ference quotient, Eq. @, which equals the slope of the in- 
clined line (green). 



Superlattices with static NDC, similarly as Gunn 
diodes, are unstable against spatial space-charge fluctua- 
tions, which result in formation of high-field electric do- 
mains. The possibility of space-charge instability in the 
superlattice Bloch oscillator was first pointed out by Kti- 
torov et al. themselves Later, it has been established 
within a simplified model by Biittiker and Thomas [l3| 
and in most consistent form by Ignatov and Shashkin 
0. 

The electric domains are believed to be destructive for 
high-frequency gain in superlattices. Nevertheless, in sc- 
ries of very interesting experiments French groups did 
observe gain at reflection of microwaves from long super- 
lattices with n ~ 10 16 cm -3 in a wide range of frequencies 
up to several hundreds of GHz [15( . It seems that physical 
reasons for these very interesting experimental observa- 
tions are still not well understood in the semiconductor 
community. 

On the other hand, for THz waves, Santa-Barbara 
group recently reported a decrease of absorption, but still 
not gain, in an array of short superlattices fl6l ] . Experi- 
mental techniques for monitoring the domains in super- 
lattices are under development |17| . 

Now we turn to consider the regime of large probe field, 
where high-frequency gain is not necessarily connected to 
the presence of static NDC. 



III. SUPERLATTICE BALANCE EQUATIONS 

Until now we have discussed electron transport in su- 
perlattices assuming the existence of a single relaxation 
time. However, experiments demonstrate that models 
employing two relaxation times are more realistic. By us- 
ing the tight-binding approximation and the Boltzmann 
equation with two scattering times, the following balance 
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FIG. 2: (Color online) Real (red) and imaginary (blue) parts 
of conductivity for (a) Edc = 0.8E cr and (b) Edc = 3E cr - 
Note that for Edc = 3E cr gain is possible almost up to lu = 
lob = 3r _1 . 



equations can be derived |8| 

d . . d 2 e 2 E(t) , . 



dt 



e(t) ~ E(t)j(t) 



-Jvj(t) 

-7e[e(*) 



-eq\ ■ 



(7) 



Here j e and j v — j e + j e i are the phenomenological scat- 
tering constants for electron energy and miniband elec- 
tron velocity respectively, 7 e ; is the scattering constant 
describing elastic scattering events, j(t) is the current 
density and e(t) is the total miniband energy density of 
electrons. Electrons are staying at the bottom of mini- 
band for e = — nA/2 and the upper edge is reached if 
e = +n A/2, where n is the density of electrons in the 
first miniband. Average electron energy in thermal equi- 
librium e eq depends on the temperature, superlattice pa- 
rameters and carrier density. 

The superlattice balance equations have been first in- 
troduced by Ignatov and Romanov [Tsj . Importantly, 
equations Q can describe both transient (t < 7~ x ) and 
stationary (t 3> 7^ ) dynamics of miniband electrons 
under the action of strong electric fields (The electric 
field is turned on at t = 0.). Earlier these equations 
have been employed in the studies of decaying coher- 
ent Bloch oscillations 0] , absolute negative conductance 
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|2Cj and efficiency of THz Bloch oscillator .21] , as well as 
in the analysis of such strongly nonlinear phenomena like 
symmetry-breaking and chaos ■ 

The solution of equations Q for static field is a scaled 
Esaki-Tsu characteristic. Its critical field is 



h/(edT eS ), 



(8) 



where r c g = (^/^WTe ) • This new definition of critical 
field is used in the following and it reduces to the earlier 
definition, if j v = 7 e = 1/r. 

In addition to parameter r c g, a parameter v — j £ /"f v 
is needed to describe the system. If v = 1, there is no 
elastic scattering i.e. j e i = 0. 

Assuming strong electric field of the form 0} , the equa- 
tions Q) have time-dependent stationary solution in an- 
alytic form if 7„ = j e = 1 Jt. In this case we find for the 
Fourier harmonics of the current I(t) = j(t)A 

IS = ^jfWET^Eded + lhu), 

I 

I h' C ° S = ^2Ji( a )[ J i+h(a) 
i 

+Ji-h(a)]l E T(eE dc d + Ihu), 
^r 31 " = X! Ji{ a )[Ji+h{a) - Ji-h(a)]K(eE dc d + Ihu), 

(9) 

where J n {x) are the Bessel functions, the summation is 
from — oo to +oo, a = eE ac d/ (hoj), Iet is given by equa- 
tion J3J) and 



K{eEd) = 7 peak - 



(eEd/ry 



The solutions of the balance equations © are the same 
as the corresponding expressions for current components 
found with the help of an exact formal solution of Boltz- 
mann transport equation with a single relaxation time 
[23|. Moreover, it is easy to see that for small ac field 
E ac (i.e. for a -C 1), /" ,cos becomes the quantum deriva- 
tive defined by eq. ©■ 

Finally, it is useful to find the expression for the dc 
differential conductivity 



dE d 



= ^jf{a)a Q {eE dc d + lhuj) 



21, 



peak 



E c 



£ J *»T7 



{eE dc d + mw) 2 /r 2 



[l + {eE dc d + lhLof/T^\ 



(10) 



where <7o (eE dc d) is the dc differential conductivity taken 
from the Esaki-Tsu characteristic. 

Formulas @ and (|TU)l will be used in the next section 
in the stability analysis of large-amplitude oscillations in 
the Bloch oscillator for the case j v = j E . They can also 
be used to check the numerical solution of the balance 
equations in the limiting case v = 1. 



Note that since in the case of static electric field the 
solution of equations was the scaled Esaki-Tsu VI 
characteristic for all u, different scattering rates do not 
give any qualitative changes to low frequency behavior. 
This is, however, not necessarily the case for higher fre- 
quencies. The quantum derivative result is valid only 
if v = 1, although it works as a good approximation, 
when v sa 1: Qualitative changes in the high frequency 
behavior can occur even for small signals, when v is sig- 
nificantly smaller than 1. The changes in large amplitude 
oscillation dynamics are even more interesting, because 
the calculations for v — 1 indicate possibility for gain 
and domain suppression ;9] . 



IV. LARGE-SIGNAL GAIN AND STABILITY 

We solved the balance equations J7J) numerically, found 
the stationary time-dependent current , and then cal- 
culated the dc current component Iff , the absorption 
jrw,cos anc j k e c[ c differential conductivity 

a uj _ aI 

dE dc 

For v — 1 these functions were found using analytic for- 
mulas (jHJl and (fTTT|> . In this case, we found excellent agree- 
ment between these two approaches, as well as with the 
results of reference 

According to the formula Q , the dc current can be cal- 
culated with the help of photon replicas of the Esaki-Tsu 
characteristic, and when the amplitude of the ac field in- 
creases, the so-called Shapiro-like steps occur in VI char- 
acteristic [Fig. |3 (v =1)]- Within these step-like struc- 
tures, regions of positive differential conductivity (PDC) 
exists [Fig. (b)]. When v decreases, the structures be- 
come less pronounced but the region of PDC becomes 
even larger. 

The PDC should be considered as one of the condi- 
tions for electric stability of the system 9] . This is a sort 
of extension of the Limited Space Charge Accumulation 
(LSA) mode of Copeland [24[, well-known in physics of 
Gunn diodes, to the case of superlattices and THz fre- 
quencies. 

On the other hand, even for a large enough ac field 
E ac the absorption can stay negative for lu < lob , as can 
be expected from the electron bunching mechanism |l0j . 
Figure 0] shows the absorption If os and the dc differen- 
tial conductivity as functions of E ac , when uit c s — 10 and 
^BT c g = 10.8. As evident from this figure, there exists a 
well-defined range of amplitudes E ac for which the posi- 
tive differential conductivity and the negative absorption 
occur simultaneously. When v decreases the range of 
amplitudes E ac for negative absorption shrinks and the 
range of amplitudes for PDC expands. 

The ranges of dc bias E dc and ac amplitude E ac , sup- 
porting PDC and gain at the first Shapiro-like step, are 
presented in Fig. [S] for three different values of v, when 
WTeff = 10. In this parameter plane, the area of gain and 
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FIG. 3: Typical current-field characteristic (wr e g = 10, 
E ac = 10-Ecr) demonstrating local structures (a) with positive 
differential conductivity (b) for u>T c g — 1 < ujbt c s < WT e g + 1. 
(Note that wsr e ff = Edc/E cr .) The local structures become 
weaker, when the ratio of scattering times decreases, but PDC 
can exist even in larger range of parameters. 



simultaneous electric stability changes with a decrease of 
v (cf. (a) and (c)). The change is very small as long as 
v > 0.4. 



V. CHARACTERISTIC TIME SCALES 

Until now we have considered only stationary trans- 
port properties. Here we briefly examine and compare 
the time scales determining the evolution of the fields 
in a superlattice. There are the following characteristic 
times. Stationary VI characteristic is established with 
the characteristic time r. (Single scattering time is as- 
sumed for simplicity.) AC field defines the time scale 
2ir/ui. Field inside an ideal (very high-Q) resonator is 




b) 



(dl /dE do )/(l peak /E cr ) 




FIG. 4: Local PDC and gain can exist simultaneously in a 
wide range of amplitudes E ac . Dependencies of J"' cos ( a ) and 
<tq (b) on probe amplitude E ac are shown, when u>T c g = 10 
and cJsTcff = cures + 0.8. If v = 1, dc differential conductivity 
is positive for 5E cr < E ac < 34E cr and absorption is negative 
for < E ac < 27E cr . This means that gain in conditions of 
electric stability exists for 5 < E ac < 27E cr . If v is decreased, 
the local structures become weaker and the magnitude of the 
absorption decreases: The gain in conditions of electric sta- 
bility occurs in a smaller range of E ac values. The changes 
are small, if v > 0.5. 



growing with the characteristic time 

= 2ee 
Tg Re[a(uj)Y 

The characteristic time for domain formation is |2f 



Td 



(11) 



(12) 



where Co is the dc differential conductivity. 

We are working with a high-frequency fields satisfy- 
ing ujt > 1. The condition of LSA, additionally to the 
requirement for a positive slope of time-averaged VI char- 
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FIG. 5: (Color online) Regions of a% > (grey) and J"' cos < 
(green) for uit c h = 10 with different values of v: (a) v = 0.1, 
(b) v = 0.4 and (c) v — 1.0. Overlapping of these regions 
(violet) gives gain without instabilities. 



acteristic (PDC), includes also the inequality 

LOT d > 1. 



(13) 



It indicates a small (limited) charge accumulation during 
every THz cycle. Condition (|13l) can be easily satisfied for 
typical superlattices. Only in heavy doped superlattices 
with very wide minibands, Td can be comparable to r 
p^ . and therefore uiTd — uit^/t) can approach unity. 

On the other hand, increasing the carrier density we 
are not only increasing the gain but also reducing the 
time scale for the growth of the field in the resonator 
(|llfl . Moreover if r g <C Td, then the resonator- mode can 
reach the minimum amplitude required for switching to 
local PDC before a domain would be developed. Such a 
scenario can potentially solve the "device turn-on prob- 
lem" for Bloch oscillator. The ratio t^/t^ depends 
only on oj and lob, while material parameters are not 
crucial. However, according to our calculations r g and 
Td are close to each other for the most interesting case 
Edc 3> E cr . This means that the growing space charge 
domains and the growing ac field cannot be handled sep- 
arately. 



VI. DISCUSSION AND CONCLUSION 

In summary, we have shown that THz gain in the 
conditions of suppressed electric domains is possible in 
fairly large regions of parameter space, which could al- 
low to build devices which can be used for generation and 
amplification of THz radiation. Large-signal gain with 
suppressed domains is preserved with an introduction of 
two different relaxation times for the electron velocity t v 
and energy r e . In particular, we demonstrated that un- 
til r v /r e > 0.5, the difference from the results obtained 
using a single r is negligible. Note that according to 
the experiments [27j t v /t e > 0.5 is a good assumption, 
although it is not valid for all superlattices (28|. Quan- 
titatively, the magnitudes of all current components and 
thus also the gain always decrease with decreasing ratio 
t v /t £ . 

There remains, however, several important problems, 
which have not been considered in the present work. The 
first one is the possible influence of boundary conditions, 
which according to the computational results of Rieder 
|29| may turn out to be crucial from the point of view of 
extended LSA mode. Hopefully, it is still possible to con- 
trol these boundary conditions at least to some extent in 
which case the LSA mode can be made to work. Second, 
we should mention that existence of LSA regime in super- 
lattices at THz frequencies is not completely proven even 
theoretically. The existence of LSA operational mode is 
rather well established both theoretically and experimen- 
tally in Gunn diodes at microwave frequencies. However, 
still no experiments devoted to the LSA mode in THz 
range are performed. From the theory side the consis- 
tent derivation of LSA conditions has been done only 
within quasistatic approximation and thus it is not valid 
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at THz frequencies. Therefore, for those who like com- 
plete proofs, the existence of extended LSA regime at 
THz frequencies still continues to be interesting theoret- 
ical problem to study. Theoretical research in these and 
related directions is in progress in Oulu. Of course, ex- 
periments would also be very helpful in order to solve the 
remaining important problems. 

Finally, in this work we have mainly focused on large- 
signal THz gain in superlattices with suppressed space- 
charge instability. A very important problem that still 
remains is how to get small-signal gain without domains. 
One possible solution is to make use of microwave pump 
to get gain at frequency multiplication in superlattices 

More full presentation of our research will be published 



elsewhere. 
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